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We investigate the effect of radiation reaction on the motion of a wave packet of a charged scalar 
particle linearly accelerated in quantum electrodynamics. We give the details of the calculations for 
the case where the particle is accelerated by a static potential that were outlined in Phys.Rev. D 70 
(2004) 081701 (R) and present similar results in the case of a time-dependent but space-independent 
potential. In particular, we calculate the expectation value of the position of the charged particle 
after the acceleration, to first order in the fine structure constant in the fi — > limit, and find that 
the change in the expectation value of the position (the position shift) due to radiation reaction 
agrees exactly with the result obtained using the Lorentz-Dirac force in classical electrodynamics 
for both potentials. We also point out that the one-loop correction to the potential may contribute 
to the position change in this limit. 
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I. INTRODUCTION 

A charged particle radiates when it is accelerated. This 
simple discovery has been the basis of a multitude of 
technological advances since. Most work, however, has 
concentrated on the radiation itself. The effect of the 
radiation, which carries energy and momentum, on the 
particle itself is very small but still present. The overall 
effect is a change in the equation of motion for the parti- 
cle. The nature of this modified equation has been sub- 
ject to much controversy since the initial work of Abra- 
ham and Lorentz. The standard equation for the change 
in the energy and momentum of the particle in classi- 
cal electrodynamics is the relativistic generalization by 
Dirac of the work of Abraham and Lorentz and called 
the Lorentz-Dirac (or Abraham-Lorentz-Dirac) equation 
[ESQ. (See, e.g., Ref. [| for a modern review.) Thus, 
if a charge e with mass m is accelerated by an external 4- 
force Fg'^j, then its coordinates x^{t) at the proper time 
T obey the following equation: 



where the Lorentz-Dirac 4-force -F^j) is given by 



^LD — 



2a, 



d^x^" 

dT^ 



dxf" 
~d^ 



dT^ dT^ 
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(1.2) 



We have let c = 1 and defined ac = e^/47r. Our metric 
is = diag (-1-1,— 1,-1,— 1). This equation of motion 
is fundamentally different in form from those which one 
usually encounters in mechanics as it is third-order. This 
implies that a third initial condition is needed in addition 
to the position and velocity. Herein lies the cause of much 
of the debate as to the physical correctness of this theory. 
It should be noted that there are many derivations of the 
Lorentz-Dirac force using a number of different methods 
besides the original work of Dirac 0] (see, e.g., Ref. 0). 




FIG. 1: Preacceleration of a charged particle 



However, many unphysical and problematic effects occur 
in the implementation. One example is the existence of 
run-away solutions for which a particle will continue to 
accelerate under its own radiation reaction. This is cer- 
tainly an unphysical effect and thus a problem for the 
theory. Alternative forms of implementing the Lorentz- 
Dirac force, including the integro-differential form, solve 
the problem of run-aways but create another problem in 
the form of pre-acceleration. Thus, in this implementa- 
tion of the Lorentz-Dirac force the particle accelerates 
before the external force causing the acceleration is ap- 
plied; this is an acausal effect which should not be present 
in classical mechanics. This effect is represented in Fig. 
n which shows the acceleration of the particle, a, against 
time, t, where a constant force is applied after time tg, 
i.e. in the shaded region. 

Additional conditions and procedures are therefore 
necessary in order to treat the Lorentz-Dirac theory as a 
a normal causal classical theory, and there arc such proce- 
dures (see, for example, the review 0] or the discussions 
in and 0)- One such procedure is the reduction of 
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order '8|, in which the Lorentz-Dirac force is treated as 
a perturbation order by order. Despite this, the debate 
over the controversial aspects of Eq. H1.2|l has continued, 
its impact muted for a couple of reasons. Firstly, the 
effect due to the third-derivative term, often called the 
Schott term, is tiny for most physical situations; in fact 
the time scale for the pre-acceleration is too small for 
any known classical interactions for it to have a measur- 
able effect. Secondly, the advent of quantum electrody- 
namics (QED), which is the more fundamental theory of 
electrodynamics, has rendered the problems concerning 
the classical Lorentz-Dirac force less urgent. The second 
point, however, naturally leads to the following question: 
"If the perturbative treatment of the Lorentz-Dirac force 
is causal and satisfactory, can it be derived in QED in 
the h limit in perturbation theory?" This question 
is the main motivation of this paper. 

In fact it was found in Refs. Hill (after an initial 
claim to the contrary) that the change of the position, 
called the position shift, due to the Lorentz-Dirac force of 
a charged particle linearly accelerated by a static external 
potential is reproduced by the h limit of one-photon 
emission process in QED in the non-relativistic approxi- 
mation. (See, e.g., Refs. for other approaches 
to arrive at the Lorentz-Dirac theory from QED.) Here 
we present in detail the generalization of this work to a 
fully relativistic particle, which was outlined in At 
the end of the paper we shall also show similar results 
with a time-dependent but space-independent potential. 
The paper is organized as follows. In Sees. ^ and IIIII 
we outline the classical and quantum models that are 
used in the paper. In Sec. IIVI we calculate expressions 
for the position expectation value for both the hypothet- 
ical non-radiating particle fSubsection IIV A|) and the ra- 
diating particle fSubsection IIV B|) before identifying the 
expression for the position shift in terms of the emission 
amplitude. This amplitude is calculated in Sec. The 
result is used in Sec. IVII to show that the energy emitted 
in the h ^ limit in QED is given by the classical Lar- 
mor formula and to compute the quantum position shift 
in Sec. IVIII The comparison between the classical and 
quantum position shifts is presented in detail in Sec. lVIIll 
and the time-dependent case is presented in Sec. llXI be- 
fore we conclude the paper in Sec. 00 



II. CLASSICAL MODEL 

We start by describing the model to be investigated in 
this paper. Consider a charged particle with charge e and 
mass m moving in one space dimension under a poten- 
tial. Let this motion be in the positive z-direction and let 
the potential be a function of the spatial coordinate z, viz 
V — V{z). We wish to analyze the change in the position 
of the particle due to radiation reaction. Thus we con- 
sider a model in which there has been a period of accel- 
eration, i.e. non-constant potential, at some time in the 
particle's history. We assume that V{z) = Vq = const. 



V{z) = Vq 



Acceleration 



V{z) = 



V 



FIG. 2: The potential V{z) and period of acceleration 



for z < ~Z\ and V{z) — for —Z-i < z for some Z\ and 
Z2, both positive constants with Z\> Zi. Thus, there is 
non-zero acceleration only in the interval (—^1,-^2) as 
represented in Fig. [21 for the case with Vq > 0. 

The external 4-force in Eq. Hl.l|l representing the elec- 
tric force corresponding to this potential is given as 



^cxt 



V'{z)dz/dT, 
V'{z)dt/dT, 

0. 



(2.1) 



The Lorentz-Dirac 4-force can be given in the following 
form: 



LD 



LD 



Fld dz/dr , 
Fld dt/dr , 

— n 



with 



Flo ^-7^(7^.). 



(2.2) 



(2.3) 



A dot indicates the derivative with respect to t. We have 
defined 7 = (1 — z^)~^/^ as usual. 

Suppose that this particle would be at z = at time 
t = if the Lorentz-Dirac force was absent. The true po- 
sition ai t = 0, which we denote 6z and call the position 
shift, can readily be found to lowest non-trivial order in 
Fld by treating the Lorentz-Dirac force as perturbation. 
The calculation can be facilitated by using the fact that 
the total energy, mdt/dT + V{z), changes by the amount 
of work done by the Lorentz-Dirac force. Thus we find 



Fi^j,{t')z{t')dt' 



d 



dz 



=5z + V'{z)5z 



where we have used 
d 



dt 



Q - 2 

mj z 



mj^z 



dt 



-V'{z) 



(2.4) 



(2.5) 
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to zeroth order, i.e. in the absence of the Lorentz-Dirac 
force. Rearranging and integrating, and then interchang- 
ing the order of integration, we obtain the position shift 
as 



1 



:dt' 



^ dz , 
dt 



(2.6) 

where vq = z\t=o is the final velocity in the absence of 
radiation reaction. We shall see that this result is repro- 
duced in QED in the ft ^ limit. 

It is useful for later purposes to find the change in the 
position shift (|2.t)|) caused by letting the particle be at 
z — ZQj^Oait — with the same final velocity vq 
in the absence of the Lorentz-Dirac force. (We assume 
that zq > —Z2 ) The time this particle spends after the 
acceleration and before getting to t = (when the posi- 
tion shift is "measured") is lengthened by zo/wo- Hence 
this position shift is obtained by using the same trajec- 
tory but defining it at = zq/vq. Hence, writing this 
position shift as 



(5zlD + 6Zc 



(2.7) 



we find 



Szc 




1 , , \ „ dz , 

^rdt' FLD-rdt. 

m J^^ yJt, j^t') [i{t')f J dt 

(2.8) 

Noting that z{t) — vq for values of t between and to we 
obtain 



5za 



,3„,2 



(2.9) 



with 7o = (1 — Vq) where -Eem is the energy emitted 
as radiation given by 



f" „ dz 

J — 00 





dt 



(2.10) 



which is the relativistic Larmor formula for one- 
dimensional motion. We shall see that the extra con- 
tribution (|2.9() to the position shift is also reproduced by 
QED in the ft ^ hmit. 



and charge e, the field is the electromagnetic field, and 
the function is the external potential which accelerates 
the charged scalar particle and is given as = V{z)Sf_io 
for the static potential. The non-interacting quantum 
electromagnetic field Af^{x) is expanded as 



d^k 

2fc(27r) 



K(k)e 



-ik-x 



4(k)e 



ik-x~\ 



(3.2) 



where k = ||k|| and the annihilation and creation opera- 
tors, a^(k) and aj((k), respectively, for the photons with 
momenta ftk in the Feynman gauge satisfy 



(k),4(k') 



-g^^(27r)32ftfc(5(k-k'). 



(3.3) 



The Fourier expansion of the non-interacting quan- 
tum scalar field is not as straightforward because its field 
equation involves the potential V{z) as follows: 



{h^[idt-v{z)f 



ftV^ 



(3.4) 



Positive-frequency mode functions can be found in the 
form 



<i>p(t, x) = 0p(z)exp 



-(p± -xx - Pat) 



Po > 0, 



(3.5) 

where x^ = {x,y) and pj^ — {p^,py). We define the 
z-component of the momentum p = p^ by the relation 
Pq — p"^ + p\ + m"^- We let p be positive (negative) if 
the mode function corresponds to a wave coming in from 
z = —00 (z = +00). The function 4>p{z) with p_L = 
satisfies 



'^^dz2 '^p(^) = [i^pi^)? ^p(^) ' 



where 



Kp{z) = y/[pQ - y(z)]2 - TO^ 



(3.6) 



(3.7) 



If p^ 7^ 0, then the mass term m? has to be replaced by 
rri^ + p\. However, we need only the cases with pj^ = 
and p > in our calculations. The WKB approximation 
for the function 4>p{z) with p > 0, which will be useful 
later, reads 



(j)p{z) 



Kp{z) 



exp 



Kp(C) dC 



(3.8) 



III. THE QED MODEL 

We now turn our attention to the quantum field theo- 
retic model. The corresponding Lagrangian density is 



C 



1 



1 



(3.1) 



where F^j,^, = d^A^, - d^^A^ and Dp^ = iV^/ft. The 
field (/5 describes a charged scalar particle with mass m 



where we have required (/'p(O) = 1. 

The non-interacting quantum charged scalar field is 
expanded in terms of these mode functions as 



</j(.t) = ft 



d3p 



2po(27rft)3 



[A(p)$p(.x)+i?t(p)cl>;(x)] 



(3.9) 

(This expansion needs to be modified if there are modes 
which are totally refiected by the potential — this is the 
case if Vo 7^ — or if there are modes which fail to 
reach \z\ = 00. However, such modifications will not 
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be relevant since the affected modes will not enter our 
calculations.) The operators A{p) and B{p) satisfy 

[^(p),^t(p')] ^ [B(p),Bt(p')] ^ {2nhf2pod{p-p'), 

(3.10) 

if one normalizes the mode functions appropriately, with 
all other commutators vanishing. The WKB mode func- 
tions $p(t,x) with 4>p{z) given by Eq. (|3.8|l are correctly 
normalized since it satisfies the usual normalization con- 
dition |$p(i,x)| = 1 in the region (— Z2,-t-oo) where 
Viz) = 0. 

Now, in the interaction picture, if there is one scalar 
particle with momentum p in the initial state, then the 
final state to lowest non-trivial order in e is a superposi- 
tion of a state proportional to the initial state and states 
with one scalar particle and one photon. Thus, the state 
A'f(p)|0) evolves as follows: 

A\p)\Q) ^ [1 + iT{p)]AHp)\0) 

i r (fk 

'^hj 2fc(27r) 



3^^(p,k)at(k)At(P)|0). 

(3.11) 



The Tip) is the forward-scattering amplitude of order 
coming from the one-loop diagram, which we do not 
evaluate explicitly. The emission amplitude A^{p, k) will 
play a central role in our calculations. Note that the mo- 
mentum hk of the photon is of order h because the wave 
number k rather than the momentum has the classical 
limit. The momentum P of the charged particle in the 
final state is determined using conservation of the trans- 
verse momentum and energy: 



Px - P± 

Po = f 



hk. 



(3.12) 
(3.13) 



The z-component P = is determined by the on-shell 
condition 



P' 



We shall consider the initial state \i) given by 



d3p 



:f{v)A\pm, 



(3.14) 



(3.15) 



where the function /(p) is sharply peaked about a mo- 
mentum in the positive z-direction with width of order h. 
We also require that the p-derivative of /(p) has width 
of order fi. This wave packet corresponds to the classical 
particle considered in the previous section. We make a 
further assumption that the WKB approximation (|3.8|) 
is valid for the momenta p where the function /(p) is 
not negligibly small. The WKB approximation is known 
to be applicable if the wavelength stays approximately 
constant over many periods (see, e.g. Ref. [l5|). Since 
the "time-dependent wavelength" is 2'Kh[Kp{z)]~^ , where 
Kp{z) is defined by Eq. H3.7|l . this condition is 



d_ 

dz 



Kp(z) 



< 1 



(3.16) 



which is automatically satisfied as long as Kp(z) and ^^(z) 
can be regarded as quantities of order . 

The normalization condition (i | i) = 1 implies 



d^p 



■|/(P)I' 



(3.17) 



(27rft)3 ' 

because of the commutation relations H3.10|l . This equa- 
tion shows that the function /(p) can heuristically be 
regarded as the one-particle wave function in the mo- 
mentum representation. In Ref. TTl it was assumed that 
/(p) was real for simplicity. This assumption will be 
dropped in this paper. 



IV. POSITION EXPECTATION VALUE 

In this section we derive a formula for the position shift 
in terms of the emission amplitude A^{p,\s.). We closely 
follow Ref. 0, correcting a few errors. In particular, 
the fourth term in Eq. H4.32|l was missing in Ref. [l0|. 

A. Non-radiating particle 



Define the charge density by 



p{x) EE - : ip^dt^p - dtif^ ■ if 
n 



(4.1) 



where : . . . : denotes normal ordering. The operator p{x) 
is the i-component of the conserved current 



n 



(4.2) 



We let the potential satisfy |Vo| < 2m, thus precluding 
the possibility of scalar-particle pair creation. Then, the 
charge density p{x) coincides with the probability density 
for the particle if there is only one charged particle in the 
state. Hence the expectation value of z is 



d^xz{p{t,x)) 



(4.3) 



We shall first need to consider the expectation value of z 
for the hypothetical scalar particle that interact with the 
external potential but not with the quantum electromag- 
netic field. This will be the benchmark against which we 
can define the position shift due to the radiation. Hence, 
we must evaluate the expectation value of the probability 
density in the absence of radiation and use Eq. 1)4. 3|l to 
find (z). The initial and final states in this case are the 
same. Thus, using the commutation relations (|3.1U|) for 
the expectation value of the charge density (|4.1|l in the 
state \i) given by H3.15|l . we find 

(P(i,x)) 



2p'^{2^hf J V2p^(2^^)^" 
$;,(t,x)5t$p(t,x) - dt%.{t,^) ■ $p(t,x)] 



(4.4) 
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The i-dependence of the mode functions is given by 
$p(t, x) (X e~*Po*/''. Hence, we find 




r/*(p')/(p) 



27rft)3- 

$;,(i,x)ci>p(t,x). 



(4.5) 

We arrange the wave packet so that at time t = it is 
located far into the region where the potential vanishes. 
Then we may approximate the mode functions as follows: 



$p(0,x) 



(4.6) 



After substituting Eq. (|4.5II in the position expectation 
formula H4.3|l we use the approximation H4.6() to write 



2$p(0,x) 



'*^^*p(0,x). 



(4.7) 



We can now use integration by parts for the variable p 
in Eq. (|4.3|l to obtain the expectation value of z for the 
hypothetical non-radiating particle, which we denote by 

zq, as 

ih r (fp' f d^p 




/*(p')/(p) 



J d3x$;,(0,x)$p(0,x). 



(4.8) 



Substituting the approximation H4.6|l for <i>p gives a delta 
function: 




X {2nh)^5{p-p') . 



(4.9) 



Performing the p'-integration after explicitly working out 
the p-derivative, we find 

dp 



Zq ~ ih 



,,nfnv)^Jiv) (4.10) 



or, by integration by parts. 



ih 

Zt) = — 

" 2 



dh 



{2TThy 



■r(p) dpfip)^ 



(4.11) 



where dp=dp — dp- This result can be interpreted as 
the expectation value of the position operator ihdp in 
the p-representation of the one-particle wave function. 
An example of a function /(p) satisfying Eq. (|4.11|) is 
/(p) = /j?Xp)e~'^''^^'' with the function /i?(p) being real 
though it is not necessary to make this assumption. 



B. Radiating particle 

Now we consider the case of the radiating particle. The 
final state resulting from the initial state \i) can be found 
from Eq. H3.11|l as 



l/> = 



d^p 



+ - 



2po(27rn)3 



F(p)At(p)|0) 

rf3p 



where = pj 
have defined 



2fc(27r)3 J V2^(27r/i)3 
xG^(p,k)at(k)At(P)|0), (4.12) 

h'k± and Pq = Po ^ ^-^ as before. We 



F(p)^[l + z.F(p)]/(p), 
G''(p,k) = ^^(p,k)/(p). 



(4.13) 
(4.14) 



One can heuristically regard the function F{p) as the 
one-particle wave function in the zero-photon sector in 
the p-representation and the function G'''(p, k) as that 
in the one-photon sector with a photon with momentum 
hh in the P-representation. By introducing the definition 



d3p 



/2^(27rn)3 
the final state can be written 

d3p 



^^(p,k)/(p)At(P), (4.15) 



I/) = 



'2^(27rn)3 

d3k 



F(p)At(p)|0) 



2fc(27r) 



Ct^'(k)at(k)|0). (4.16) 



We first derive the relation between the imaginary part 
of ^(p) and the scattering probability which results from 
unitarity. Recalling the definition H4.13|l and using the 
unitarity of time evolution, we find 



l = (/l/) 



_fp_ 

(2^ 



(l-2Im.F(p))|/(p)P+Pem, 



to first order in e^, where the emission probability Vcm 
is given by 



d^k 
2fc(27r)3 

1 r d^k 

h J 2fc(27r)3 
d3p' 



(0|C^(k')Ct^(k)|0) 



d3p 



/2p^(27rfi)3 
x(0|A(P)At(P')|0). 



%/2R^(27r?i)3 
/(p')-4^(p',k) 



/*(p)^;(p,k) 



(4.18) 



To perform the p'-integration recall first that 

(0|A(P)At(P')|0) = {2TThf2Po6{P - P') . (4.19) 
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We need to integrate this delta function with respect to 
p' by relating it to S{p — p'). We note that 

5{P^ - PI) = 5{p^ - p1) (4.20) 

because = p^ — fik±, and 

6iP-n = ^Sip-p'). (4.21) 

By using these formulae in Eq. (|4.19|) wc find 

(0|A(P)At(P')|0) = ^^{2nnf2poS{p - p') . (4.22) 
Pq dP 

By substituting this formula in Eq. H4.18|l we obtain 

^ii^x;,(p,kM-(p'.k)5i* 



2fc(27r)3 



Po dP 



(4.23) 



Eq. (|4.17|) must hold for any function /(p), so we must 
have 

(4.24) 

The contribution of the one-photon-emission term to 
(p(t,x)) is 



^^^(0|C^(k)p(t,x)Ct(k)|0). 



(4.25) 

The integrand here is obtained from {p{t, x)) for the non- 
radiating particle if we replace the state \i) in Eq. (|3.15() 

by 



Si/2ct(k)|0)=z 
where 



/2Po{2T:h) 



3.g,(P)At(P)|0), 



g,(P)^?.-V2^,(p,k)/(p)./^^. 

\j Po dP 

Hence we conclude from Eq. (|4.11|l that 

h J d3xz(0|C^(k)p(i,x)Ct(k)|0) 

ih f d^P 



(4.26) 
(4.27) 



(27r)3 



g^*(P)apg^(P) 



(4.28) 



This implies that the one-photon-emission contribution 
to (z) is 



(^>c 



d^k 



d^P 



2 J 2fc(27r)3 J {2TThf 
X [/*(p)^''*(p,k)] Hp [/(p)^^(p,k)] 



Po / dp 



PO \dP 



(4.29) 



We convert the integration variable from P to p by using 
d'^P = {dP / dp)d^p and the P-derivative to a p-derivative 
by using dp = {dp/dP)dp to find 



dh 



2.1 (2nh)^ 
d^k 



rip) dpfip) 



Pq f dp 



2fc(27r) 
I f d^p 



3^'^*(p,k)^,(p,k)^^ 



2.1 (2nh)3 
d^k 
2k(2nf 



I/(P)P 



^^*(p,k) ap^p(p,k). (4.30) 



We have dropped the factor {Po/po){dp/dP)'^ from the 
second term for the following reason. It will turn out 
later that the emission amplitude and its p-derivative 
are both of order W . Thus the second term in the above 
expression for {z)f.m. is of order W . Since Pq = po and 
dp/dP = 1 to order hP, we can drop the factors of Po/po 
and dp/dP from the second term. However, one needs 
to keep these factors in the first term, which is of order 

h-\ 

From the expressions for the final state (|4.16|l it is clear 
that the contribution to (z) from the process without 
photon emission is obtained from Eq. (|4.11|l by replacing 
/(p) with P(p) = [1 + z.F(p)]/(p). Thus 

(^)for = -0^{2^\fip)\%Rc:F{p) 



[l-2Im.F(p)]r(p) dpfip)} (4.31) 



to first order in e^. (Recall that the forward-scattering 
amplitude ^(p) is of order e^.) We can now add the 
two contributions (z)cm and (z)^^ together to find the 
position expectation value: 



2 



d^p 
(2^ 
d^p 



(27r/i) 
I I d^p 



2.1 (27rft)3 
d^k 



/*(p) dp fip) 

^\f{p)\'dpReJ^{p) 
I/(P)P 

-3^^*(p,k) dp Apip,k) 



2k{2iT 

i/(l|^^>)^^^(P) 
d^k 

.^^*(p,k)^^(p,k) 



2fc(27r)3'' v^'-v-MV 



dp f dp ^ 
^ Po dP \dP 



(4.32) 



The unitarity relation (|4.24|) has been used to eliminate 
2ImJF(p). Each term in Eq. (|4.32|) can be interpreted 



7 



as follows. The first term is the position expectation 
value zq for the non-radiating particle. The second term 
is the contribution from what can be regarded as the 
one-loop quantum correction to the potential. Therefore, 
we identify the sum of the third and forth terms as the 
position shift to be compared with the classical position 
shift (^Zciass given by Eq. (|2.7|) . Using Eq. (|4.11|) and 
the assumption that the function /(p) is sharply peaked 
about a momentum in the positive z-direction with the 
width of order h, we find these two terms in the h ^ Q 
limit as 



EMISSION AMPLITUDE 



6zg ~ 5z„i + 5z, 



q2 



(4.33) 



where 



5zqi = -- 



0Zq2 = - 



2 J 2fc(27r)3 
xy^'^*(p,k) y^p(p,k), 

h J 2fc(27r)3 



(4.34) 



dp 



x^''*(p,k)^^(p,k) yjp-^j '(4-35) 

where the momentum at which the function /(p) is 
peaked is now denoted simply by p. We have used 
the fact that dp/dP — 1 is of order h to drop the fac- 
tor {Po/po){dp/dP) in Szq2. We shall demonstrate that 
the quantum position shift Szq in the h —^ limit is 
identical with the classical counterpart Szdass by show- 
ing SZqi = SzlB and SZq2 = SZextra- 

Let us first examine the latter equality. To this end 
we need to find an expression for dp/ dP in terms of p 
and k. The energy conservation equation pq — Pq = hk 
gives a one-to-one relation between p and P for a given 
k after letting = p^ = because these are of order 
h^. Then we find 



dp 
dP 



1 



hk. 



By using this formula in Eq. (|4.35|) we obtain 



0Zq2 — 9 ^cm ; 



p'^Po 



(4.36) 



(4.37) 



where 



d^k 



2fc(2-)3^-^'*^P'^)-^''(P''^) (^-^^^ 

is the expectation value of the energy emitted as radia- 
tion. By comparing Eq. H4.37|l with Eq. (|2.9|l . it can be 
seen that the equality 6zq2 — (Jzoxtra will hold if 

2a /■° 

£cm = ^ h'zfdt, (4.39) 

•J J-oo 

which is identical to the relativistic generalization of the 
classical Larmor formula. To show this equality (and 
also Szqi = Szi^Y)) we need to find the ?i — > limit of 
the one-photon-emission amplitude ^^(p, k) to which we 
now turn. 



The one-photon-emission part of the evolution of the 
state was represented earlier in Eq. (|3.11|) by 



At(p)|0) 



^^^-(p,k)at(k)At(P)|0). 



(5.1) 

The evolution in perturbation theory in the interaction 
picture is generated by the interaction Lagrangian den- 
sity as follows: 

AHpM -^■■■ + JiJ d^^Ci{x)A^{pm . (5.2) 

We can take the inner product of the two expressions 
(|5.1() and H5.2|l with the state (0|aj^(k)A(p') and equate 
them in order to find the amplitude ^^(p,k). This pro- 
cedure leads to 

-2hp'o{27rhf6ip' -P)AAp,'k.) 

= J d4x(0|a,(k)A(p')£/(a:)At(p)|0) . (5.3) 

Integrating over the variable p', with the appropriate 
measure, and rearranging, we find that .4p(p, k) is given 

by 



A'(p,k) 



d^p' 



h J 2p'g{2'Kh)^ 

X J d^x{0\a^{k)A{p')jri{x)A^{p)\0) . 

(5.4) 

The Ci (x) is obtained from the Lagrangian density H3.1|) 
as follows: 

Ciix) = -jA^ : [ip^D^^ip - (i?^^)V] : • (5.5) 

Recall that D^^ = 5^ + iV^, where = 6i^qV{z). By 
substituting Eq. I|5.5() in Eq. 1)5.4(1 and using the commu- 
tation relations (|3.3|) and H3.10|l we find 



A^,{p,^k) 



ieh 



d^p' 



d^xe'^-- 



2p'^{2T:hf 
$;,(x)i?^$p(x) - {D^<^^,{x))^ <i>p(x) 



(5.6) 

where the mode functions ^'p(x) are given by Eq. (|3.5|l . 
As we mentioned before, we assume that the WKB ap- 
proximation 1(3. 8|l is valid for the initial mode function. 
Since the momentum hk of the emitted photon is of or- 
der h, the WKB approximation can also be used for the 
mode functions for the final state. 

Since the transverse momentum p^ is assumed to be 
of order h, the final transverse momentum px — ^kj^ 
is also of order h. This means that we do not need to 
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consider the x- and y-components .4p(p, k) because they 
are smaller by a factor of h compared to the t- and z- 
components. Using Dq = 9t + iV{z)/h we obtain 



A(p,k) 



X 4>l,{z)4>p{z)[p^+p'^~2V{z)] 
X e*'°''^e^^['''o^^"'*^'p'^~P-^^'''-^]/'' . (5.7) 
Integration with respect to t, x and y gives 
dp' 



A(p,k) 



X [po+Po-2l^(z)],^;,(z)(/.p(z)e- 
X(5(po + hk-po) . 



(5.8) 



We convert the p'-integration to pp-integration by noting 
that dp' Ip'q — dp'^/p' and find 

A(p,k) = -e j dz^^^±P^^^rp{z)M^)e-'''\ 

(5.9) 

where Pq — Po ~ with P = \/ Pq ~ m? as before. We 
let Pq — Po a-iid P = p in Eq. I|5.9|l because the differences 
^0 ~ Po and P — p are of order h. Thus, we obtain 



A(p,k) 



-e / dz^^-^^(l)*p{z)<l)p{z)e''^'\ 



(5.10) 

Proceeding similarly for the z-component ^^^(pjk) we 
have 



A(p,k) = J dz 



,* f ^ d<t)p{z) d(f)*p{z) ■ 
^P^^^—^ ^'^p(^) 



dz 



(5.11) 



To simplify the expression for the emission amplitude Af_i 
further, we need to use the explicit form of (t>p{z) given 
by Eq. We note first 



(l>*p{z)(l>p{z) 



Pp 



Kp{z)Kp{z) 



X exp 



[Kpiri) - Kp{i])]d-q 



(5.12) 



Now, to lowest order in h we have 



Kp[z) - Kp[z) = — (po - Po) = — nk . 

opo Hp{z) 

(5.13) 

Note that Kp{z) and pq — V{z) are the z- and t- 
components of rndx^/dr, respectively, of the classical 
particle with final momentum p and vanishing transverse 



momentum. Hence, Vp{z) = Kp{z)/[po — V{z)] is the ve- 
locity of this classical particle. Thus, 



Kp{z) - Kp{z) = 



hk 

Vp{z) 



(5.14) 



By substituting this formula in Eq. H5.12|l and letting 
P = p and Kp(z) = Kp{z), because the differences P — p 
and Kp{z) — Kp{z) are of order H, we find 



(t)*p,{z)(j)p{z) 



Kp{z) 



exp 



dz 



M^) 



P ^ikt 



Kp{z) 



(5.15) 



where the time t is defined by dz/dt — Vp{z) with the 
condition t = at z = 0. Hence 



A(p,k) 



dz 



Akt — ik^ z 



Vp{z) 

^^^zkt-ik'z ^ 



It is convenient to define the variable ^ by 

^ = t — z cos 9 , 
cos6' = k^/k. 



Then 



A(p,k) 



di 



dt 
dl 



(5.16) 



(5.17) 
(5.18) 



(5.19) 



Next we turn our attention to the z-component given by 
Eq. ()5.11(l . We note that, to lowest order in h, we can let 



^-,!!4f)^^(,). (5.20) 



dz 



Hence 



A(p,k) = e j 



p+oo 

1 dz 

J —oo 


Kp + Kp 

2P ' 


r+oo 

/ dz 


^ikt—ik^ z 






r+oc 

/ 




t Kp(z) 


= Kp(z) 



(l>*p{z)(l>p{z)e' 



~ik^ z 



(5.21) 



Eq. (|5.15|) . This formula and Eq. (|5.19|) can be combined 



as 



^'^(p,k) = -e / de^e^'^-« 
J-oo d£, 



(5.22) 



where x^ is the classical trajectory with final momentum 
p which passes through (t, z) = (0,0). (The minus sign 
for the z-component is due to the raised index.) This 
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emission amplitude is identical with that for a classical 
point charge passing through (0,0) TIjI. 

The expression H5.22|l is ill-defined since dx^/d^ re- 
mains finite as ^ — > ±00. Therefore, we introduce 
a smooth cut-off function xiO which takes the value 
one while the acceleration is non-zero with the property 
linv^+^ yf£) = 0. Thus, 



+ 00 i n 



dC 

In the end we take the limit xiO 
such a way that J^^ [x'iO? - 



(5.23) 



1 at any given ^ in 



VI. DERIVATION OF THE LARMOR 
FORMULA 

In this section we derive the relativistic generaliza- 
tion of the Larmor formula for one-dimensional motion, 
Eq. H4.39|l . thus completing the demonstration of the 
equality 5zq2 = Szc^tm- It is convenient to use the fol- 
lowing form of the emission amplitude obtained by inte- 
gration by parts: 



^'^(p,k) 



le 
T 



+ 00 



d£,^ dC 



(6.1) 

where we have used the condition that x(0 = 1 if 
d^x'^/d^-^ 7^ 0. By substituting this formula in Eq. 
we obtain 



dn 



dk 

16^ 
dxi" 



de 



d^ 



(6.2) 



where dfl is the solid angle in the k-space. We extend the 
integration range for k from [0, -l-oo) to (— cxd, +00) and 
divide by two. Then using e''^^?-^') dk ^ 2n6{^-^'), 
we find 



167r2 



dn 



+ 00 



de de 



dc 



(6.3) 



The second term tends to zero in the limit x(0 ~* 1 due 



to the requirement J 
in this limit 



+00 



47r 



+00 



[x'(^)]^d^ 0. Hence, we have 
d^Xa d'^x^ 



d£_ dfl- 



(6.4) 



de de 

where ac = /An as before. Now, one can readily show 
that 

(6.5) 



d^x" 



de 



d^ d^x^' 



dt dt^ 



d^i da;'' 
dt^'df 



By substituting d^/dt = 1 — icos0 we find 



de 

dH d^z 



(1 — icos6')3 



and hence 



de 



de 



COS( 



e sin^ I 



dt 



de de (1 - i COS 61)5 d^ 



(6.6) 
(6.7) 

(6.8) 



By substituting this formula in Eq. I|6.4|l we obtain 
Eq. (|4.39|) . thus demonstrating the equality 5zq2 = 

ra- 



VII. QUANTUM POSITION SHIFT 

Our next task is to show that 5zqi = <5zld- This will 
establish that the classical position shift 6zc\ass is equal 
to the quantum one 5zq in the fi. — > limit. 

In the product y4^*(p, k)i9py4.^(p, k) in Eq. (|4.34|l we 
use Eq. (jEU for A^'* and Eq. (jF!^ for dpA,,. Proceeding 
in a way similar to that led to Eq. H6.4|l we find 



6z, 



91 



-— / dn 

An 



d^ 



dxf" 

de dp\~d( 



1 a_ / dx^da^\ d_ 2 
'Adp \dS, d£, ) dC ^^^^'^ 



.(7.1) 



Next we demonstrate that the second term, which still 
contains the cut-off function, vanishes when integrated 
over the solid angle and ^. Noting that 



dx^ dx^ 
~d^~d^ 



dx^ dXf^ 
dr dr 



(7.2) 



we find by integration by parts that the contribution of 
this term is proportional to the following integral: 



dn 



d 

dC 



_ d_ fdrT 

d'i\d~p\d~i. 



(7.3) 



Since 3/ dp is taken with ^ fixed, the ^- and p-derivatives 
commute. Hence this integral is equal to 



d_ 

dp 

d_ 

dp 



dn 



de 



dn 



dr 

df 



(7.4) 



where we have used the fact that the ^-derivative of 
(dr/de'^ is non-zero only if the acceleration is non-zero 
and, therefore, only when the cut-off function x{0 equals 
one. Since 



de) " (di/dT)2 



(dz/dr) 



= An , 



(7.5) 
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we have / = 0, and thus the contribution from the second 
term in Eq. H7.ll) vanishes. The remaining term gives the 
main contribution to the position shift due to radiation 
reaction in QED in the h ^ hmit as 



d^^ dp V 



(7.6) 



To compare Szqi given by this equation with (Jzld 
given by Eq. (|2.6|l we need to find an expression of Szqi 
in terms of t rather than ^. Thus, we need to ehminate 
the variable ^ using its definition ^ = i — z cos 9. The two 
components of the second derivative d^x^/d^^ are given 
by Eqs. 1)6. 6|l and 1)6. 7|l . Let us consider the second factor 
of the integrand of H7.6|l . i.e. the p-derivative of dx^/d^. 
By interchanging the order of integration we obtain 



d_ (dx^\ _ Md_ 



(7.7) 



where the subscript ^ indicates that the partial derivative 
with respect to p is taken with ^ fixed. Then by differ- 
entiating the equation t ~ z cos + ^ with respect to p 
with ^ fixed, we find 



dp], \dp, 



(7.8) 



Combining this equation with Eqs. H6.6|) and (|6.7() . which 
gives d^x^/d^^, we can write the integrand in Eq. H7.6(l 

as 

d^x'' d fdxi,\ z d f dz\ . o„ 

— snr6'. 



dC2 dp y 



{l-zcosOfdt \dpj ^ 



(7.9) 

We substitute this expression in the quantum position 
shift H7.6|l and change the integration variable from ^ to 
t and integrate by parts. The result is 



6Zql 



47r 



dn 



dt- 



dz 
dp 



dt 



sm 



(1 — i cos ( 



(7.10) 



where we have also changed the integration range from 
(—00,-1-00) to (—00,0] because z — Q for [0,+oo). Fi- 
nally, we need to relate (dz/dp)^ to {dz/dp)t- By re- 
garding z as a function of t and p we have 



dz\ 

dz = zdt + ( dp . 

dpJt 



(7.11) 



By substituting d^ = d^ -I- cos 9 dz in this equation and 
rearranging, we obtain 



dz 



Thus, 



1 — i cos 



1 



1 — i cos 9 \dp J f 



dz 



dp. (7.12) 



Using this expression in the position shift H7.10|l we have 
dzqi ■- 



47r 





•0 / 


dfl J 


dt 1 




-00 \ 


dz 


sin^ e 


dt (1 


— z cos 



3z 



2 sm 9 cos ( 
(1 — zcoa9 



(7.14) 



We can now perform the integration over the solid angle 
with the following result: 



SZql 



2qc 
3 



dt 



dz^ 
9P/ t 



dt 







2ac d , 



dz^ 
9P/ t 



dtF^B 



(7.15) 



where we have used the expression of i^LD given by 
Eq. H2.3I) . This formula will be used in the next section 
to show that Szqi — Jzld- We note here that Eq. (|7.15|) 
is valid for an external force of any form provided that 
the photon emission amplitude is given by Eq. (|5.22l) . 



VIII. COMPARISON OF CLASSICAL AND 
QUANTUM POSITION SHIFTS 

Comparison between Eqs. (|7.15|) and (|2.6|l shows that 
the equality Szqi = (5zld follows if 



dt' 



dpJt m''^' Jo ^Ht')[z{t')] 



■dt' . 



(8.1) 



This equation can be demonstrated as follows. Since the 
energy is conserved, we have 



^Jrrfi {dz/dr f + + V{z) = ^/p^ + m? , (8.2) 
(2 (Vp^ -t- m2 - F(z)) 



and, hence 



dz 
'dt 



-,1/2 



1 — m/ 



(8.3) 



By differentiating both sides with respect to p with t 
fixed, and noting that 



P. 



/V? 



\/ p'^ + m? — V{z) — , 



we obtain 
d { dz 



dt \dp J f m^^z 



vo - V'{z) 



dz 
dpJt. 



By substituting the formula V'{z) 
Eq. |j2!3] in Eq. (jHSJ we find 



(8.4) 
(8.5) 

. (8.6) 

-mj^z [see 



r1 



1 /f) 



1)n 
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Then by integrating this formula, remembering that z = 
Q &t t — for all p, we arrive at Eq. (|8.1|) . 

The derivation of 5zqi — (5zld above uses an explicit 
relation, Eq. (|8.3|) . between the the trajectory z{t) and 
the momentum p. However, the simple form for the quan- 
tum position shift 5zqi given by Eq. H7.15|l suggests that 
there should be a derivation which does not rely on the 
relation H8.3|l . We show next that there is indeed such a 
derivation which is valid in a more general setting. 

We consider a charged particle moving in the z- 
direction accelerated by an external force which may de- 
pend on both z and t. Thus, we consider the following 
equation of motion: 



(fz 

TO-^ = [Fo^t{t, Z) 



dt 



This equation describes a classical particle under the in- 
fluence of an external force Fcxt{t,z) and the Lorentz- 
Dirac force. The model analyzed so far in this paper is 
the special case given by i^cxt(i, z) = —V'{z). 

We first consider a solution to Eq. (|8.8|l in the absence 
of radiation reaction, i.e. with i^LD set to 0, with the 
momentum and the position being p and 0, respectively, 
at <: = 0. We define P = mdz/dT and let this solution 
be denoted 



(8.9) 



We have (z(0),P(0)) = (0,p). Let a linearized solution 
about (z, P), also with Fld = 0, be given by 



(z,P) = (z + Az,P + AP) 



(8.10) 



From the equation of motion (with Pld — 0) the quanti- 
ties Az and AP are found to satisfy the equations 



—Az = 
dt 

dt 



m-\l ~ |2)3/2^p ^ A[t)AP , 
aPcxt 



dz 



Az = B{t)Az. 



.11) 
.12) 



Let {Azs{t), APs{t)) be a set of solutions labeled by s of 
these equations satisfying {Azs{s) , APs{s)) = (0, 1). The 
linearized solutions represent a particle whose position 
at t = coincides with z(0) (see Fig. The quantity 
{dz/dp)t in Eq. H7.15|l is the rate of change in the position 
z at time t as the momentum P at t = is changed while 
the position z at i = is fixed. Thus, 



/9z 
\dp 



Azo{t) 



= Azo(i) 



Ai^o(O) 

By substituting this formula in Eq. H7.15fl we obtain 



(8.13) 




FIG. 3: The world lines for the solutions {z{t), P(t)) and 

(2(t) + A2o(t),P(t)-HAPo(t)) 



and now let (z -f 5z, P -\- SP) be the approximate solu- 
tion obtained by perturbing (z, P) to first order in Pld 
with the initial conditions {6z,6P) — (0,0) before the 
acceleration starts. The quantities 6z and 6P satisfy 



dt 
d_ 
dt 



Sz = A{t)SP, (8.15) 
6P = B{t)6z + Fi^^it) , (8.16) 



where A{t) and B{t) are defined by Eqs. (|8.11|) and 
(|8.12|) . respectively. The quantity Sz {5P) is the differ- 
ence in position (momentum) between the hypothetical 
non-radiating particle and the radiating particle. Thus, 
the quantity Sz at t — is equal to (5zld- 

The solutions to Eqs. (|8.15() and H8.16|l with the ini- 
tial conditions specified above can be expressed in terms 
of another set of solutions {Azs{t), APs{t)) of the non- 
radiating equations, Eqs. (|8.11|) and H8.12|l . as 

Sz{t) = I dsFi^u{s)Azs{t), (8.17) 

^ — OO 

SP{t) = I dsFi^u{s)AP,{t). (8.18) 

These can readily be shown to satisfy Eqs. (|8.15|) and 
(I8.1()|) by direct differentiation using the conditions 
Atz{t) = and AtP{t) — 1. The classical position shift 
due to the Lorentz-Dirac force, at time t = 0, is therefore 
given by 



Sz 



LD 



dtPLD(i)Azt(0), 



(8.19) 



SZgl = - 



dtFi,uit)Azoit). 



(8.14) 



Next a similar expression will be derived for 5zld . We 
again let (z{t), P{t)) be a solution of H8.8|l with Pld = 0, 



This expression for the classical position shift i5zld is 
similar to that for the quantum position shift Szqi in 
Eq. H8.14|) . the only difference being that the quantity 
Azo(i) is replaced by — Azt(O). In Fig. 2|we show the 
world lines corresponding to the solution (z(s), P(s)) and 
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FIG. 4: The world lines for the solutions (z(s),P(s)) and 
{z{s) + Azt(s),P{s)+APt{s)) 



the approximate solution {z{s)+Azt{s), P{s)+APt{s)) to 
the equations for the hypothetical non-radiation particle 
motion. 

Now, we only need to establish that Azo{t) — — Az((0) 
to show that the quantum and classical position shifts 
(|8.14l) and H8.19|l are equal. In fact one can show in 
general that Azs{t) = — Azt(s) for any values of s as 
follows. Given a pair of solutions {Az^-^\ AP^^^^) and 
(Az(2), AF(2)) to Eqs. and we define the 

symplectic product by 



(Az(i),AP(i)|Az(2),AF(2)) 
= Az(i)AF(2)_Ap(i)Az(2) 



.20) 



It can readily be shown that the symplectic product is 
time-independent by differentiating it with respect to t 
and using Eqs. and ^J^. With the identifica- 

tion (Az(i),AP(i)) = {Azs,APs) and (Az^^), AP^^)) = 
(Az(, APf) we obtain 



Azs{t)APtit) - APs{t)Azt{t) 
= Azsis)APt{s) - AP,{s)Azt{s) 



.21) 



1, 



Since Az^(s) = Azt{t) = and AP^s) = APt{t) 
we have Azg(t) ~ — Azj(s). Hence the equality 5zqi — 
(5zld holds. This equality and the equality 5z 
establish 5zq = Szdass [see Eqs. (I2.7|l and H4.33|l ]. Thus, 
the position shift due to QED to order in the limit 
S — > is equal to the position shift due to the Lorentz- 
Dirac force for linear motion if the acceleration is caused 
by a static potential. 

We have shown that the classical position shift (5zld 
equals 5zqi in Eq. H8.14|l for any external force which 
depends on t and z. We have also shown that Eq. (|8.14|) 
is valid provided that the photon emission amplitude is 
given by Eq. 1^5. 2^) . Thus, the equality 5zqi = i5zld 
holds if the amplitude ^^(p,k) is given by Eq. H5.22|l . 
We have demonstrated Eq. H5.22|l for the case where the 



external force fext is i- independent. The extension to 
the case where the external force depends on both t and 
z does not appear to be straightforward. However, our 
method can be applied to the case with an external force 
that depends only on t because the WKB approximation 
is very similar to that for the i-independent case as we 
demonstrate in the next section. 



IX. TIME-DEPENDENT POTENTIAL 

Let the potential be given by Vt = 0, = — ^(0 
in the QED Lagrangian density %\.\\ . We show below 
that the emission amplitude in this case is also given by 

Eq. inss- 

Since the calculation is quite similar to the time- 
independent case, we present only a brief description. 
The WKB approximation to the mode functions is given 
by 



$p(t,x) « 0p(t)exp 



Po 

C7p{t) 



exp 



(p_L • +p^z) 



apiOdC 



^ Jo 



where 



, (9.1) 
, (9.2) 

(9.3) 



if p_L = 0. (The mass term needs to be replaced by 
TO^ -I- p5_ if p± 7^ as before.) The function ap{t) is the 
t-component of the kinetic energy-momentum mdx^ /dr 
of the non-radiating classical particle at t. Note that 
the momentum p including the z-component — p is 
conserved, but the energy is not. The calculation of the 
emission amplitude follows the same pattern as before 
with the roles of t and z (and those oip and po) reversed. 

Let yl^(p, k) be, again, the emission amplitude for the 
process in which the charged particle with momentum 
p emits a photon with momentum /ik so that the final 
momentum of the charged particle is given by P = p— ?ik. 
We again assume that the transverse momenta are of 
order h and can be neglected. Then the z-component to 
order [to be compared with Eq. H5.10|l ] is 



A. 



dt 



p V{t) ~ 2 

-(j)p(t)(t>p{t) e 



Po 



(9.4) 



For the i-component [to be compared with Eq. (|5.11|) ] we 
obtain 

ieh 



At 



2po 



dt 



*pdt(l)p 



5) 



The product of mode functions in Eq. (|9.4|) is 



^*p{t)^p{t) 



PoPo 



ap{t)ap{t) 



X exp 



(ap(C)-ap(C))dC 



(9.6) 



13 



To lowest order in h we can change Pq to po and ap{t) to 
ap{t) in the pre- factor as the differences are of order h. 
Recalhng the classical momentum conservation equation, 
mdz/dT + V{t) = p, and the formula rndt/dr = ap{t), 
where z{t) is the classical trajectory with z(0) = 0, we 
obtain to lowest order in h 



ap{t) - (7p{t) 



dp 



(p-P) 



dz , 
'dt' 



(9.7) 



Thus the product of the mode functions in question can 
be written 



Po 



4'*p{t)4)p{t) = — exp 



* , dz . 
dt—hk^ 



dt 



Po 



= — exp {—ik^ z) 



(9.8) 



This gives the z-component of the emission amplitude as 



A(p,k) 



dt 



di 



P-V{t) 

CFpit) 

dz 



ikt — ik^ z 



(9.9) 



where ^ = t — zcosO and cos 6* = fc^/fc as before. To 
calculate the t-component of the amplitude, note first 
that to lowest order in h we have 

Hp,{t)dt^p{t) = -iap{t)$*{t)^p{t) 

= —ipocxp{—ik^z), (9.10) 

where Eq. H9.8|l has been used. By using this equation in 
Eq. H9.5|l we obtain 



A(p,k) 



dt 



Akt — ik^ z 



di%^^^. (9.11) 
di 



Eqs. and can be written yt^(p,k) = 

—e J d^ {dx^ I dS^e^^^ . As we mentioned before this equal- 
ity is sufficient to conclude that 5zqi — (5zld- 

The equality 5zq2 = ^-^oxtrai which is necessary to es- 
tablish the equality of the classical and quantum position 
shift, is trivial in this case. Since the system is transla- 
tionally invariant, the classical position shift is invari- 
ant under the shift in the position, (z) = ^ zq, of 
the non-radiation hypothetical particle at t = 0. Hence 
(^-^cxtra — 0. Furthermore, since P ^ p — hk^, we have 
dp/dP = 1 and hence Szq2 = in Eq. H4.35|l . Thus, we 



have Szq2 = Sz^ 



0. Hence, Szq 



6z, 



class: l-B. 



the 



position shift from the Lorentz-Dirac force agrees with 
that in QED in the limit h —^ to order for a space- 
independent potential as well. 

The equality Szqi = Jzld can be verified by a more 
explicit calculation as in the case with a static potential. 
Note that the momentum conservation equation reads 



dt 



[m-fz + V{t)] = i^LD . 



(9.12) 



Hence 



S{mjz) = mj''j^iSz)^ J F^^{t')dt'. (9.13) 
Thus, the classical position shift is 

(JZclass = <52LD - - y" ^'^^') ^LDdi , (9.14) 

where we have interchanged the order of integration. 
Thus, what we need to establish is 



1 



dz 



dpJt Jo 



dt 



(9.15) 



for the hypothetical non-radiating particle. The momen- 
tum conservation for this particle in the z-direction reads 



dz 

m— + V(t)^p. (9.16) 
dr 

Hence, with the condition z = at t = 0, wc find 



I [p-vitw 



-1/2 



dt. 



(9.17) 



By differentiating this expression with respect to p and 
using p — V(t) — mdz/dr and y/[p — V{t)]'^ + w? = 
rndt/dr, we indeed obtain Eq. H9.15|l . 

X. CONCLUSION 

In this paper we showed that the change in position 
due to the radiation reaction of a particle, which we call 
the position shift, according to the Lorentz-Dirac theory 
in classical electrodynamics is reproduced by the h ^ 
limit of QED with a scalar charged particle. The calcu- 
lation was performed to lowest non-trivial order in for 
wave packets linearly accelerated for a finite time by ei- 
ther a purely space-dependent or time-dependent poten- 
tial. The quantum wave functions constituting the wave 
packets were approximated using the WKB functions for 
each potential. The agreement between QED and the 
Lorentz-Dirac theory was demonstrated by concrete cal- 
culations which gave the position shifts explicitly. We 
also showed by a general argument that this agreement 
holds as long as the emission amplitude coincides in the 
h ^ limit with that for a classical point charge. 

This work extended the corresponding work in the non- 
relativistic approximation [olllOjl and provided the details 
omitted in Ref. Jll| as well as the analysis for the time- 
dependent but space-independent potential. It will be 
interesting to generalize our results to the case with a 
particle moving in three space dimensions. It is also im- 
portant to estimate the one-loop correction to the poten- 
tial. In the (unrealistic) model in which the acceleration 
is achieved by a space-dependent mass term fl^ . this 
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correction is of order h~^, thus overwhelming the contri- 
bution from the Lorentz-Dirac force. If this correction 
turns out to be of order in the models studied here, 
which are more realistic, then the one-loop correction 
will be as important as the Lorentz-Dirac force in deter- 
mining the motion of the charged particle. It will also 



be interesting to see whether our results can be justified 
when more than one photons are emitted: our results, 
as they stand, are logically consistent only if the emis- 
sion probability Vcm given by Eq. 14.23|1 is much smaller 
than one so that the process is dominated by onc-photon 
emission. These issues are currently under investigation. 
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